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Abstract 

Hamiltonian integration methods for the Vlasov-Maxwell eqnations are developed by a 
Hamiltonian splitting techniqne. The Hamiltonian fnnctional is split into five parts, i.e., the 
electrical energy, the magnetic energy, and the kinetic energy in three Cartesian components. 
Each of the snbsystems is a Hamiltonian system with respect to the Morrison-Marsden- 
Weinstein Poisson bracket and can be solved exactly. Compositions of the exact solntions 
yield Poisson strnctnre preserving, or Hamiltonian, integration methods for the Vlasov- 
Maxwell eqnations, which have snperior long-term fidelity and accnracy. 
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The dynamics of charged particles in a plasma interacting with the self-consistent elec¬ 
tromagnetic fields can be described by the Vlasov-Maxwell (VM) equations. In modern 
plasma physics and accelerator physics, numerical integration of the Vlasov-Maxwell equa¬ 
tions is an important tool for theoretical studies, and varieties of numerical algorithms have 
been developed. Recently, geometric integration methods which are designed in the 

spirit of preserving the intrinsic structures of a dynamical system, have been developed for 
plasma physics applications js-l^. By preserving properties such as the Poisson structure 
of a Hamiltonian system and the invariant volume form of a source-free system, geometric 
integration methods usually generate numerical results with superior long-term behavior 
compared to other methods jj, l2^, and are thus more suitable for large-scale, long-term 
simulations. It is known that the Vlasov-Maxwell system is a Hamiltonian system with 


respect to a Poisson bracket 


21 


23|. In a recent paper 18|, Crouseilles, Einkemmer, and 


Faou proposed an innovative Hamiltonian split ting method for the Vlasov-Maxwell equa¬ 


tions based on a bracket hrst suggested in Ref. 


21| . This splitting scheme results in three 


solvable subsystems, whose exact solutions can be combined into algorithms for the Vlasov- 
Maxwell equations that preserve the structure specified by the bracket. However, it has 


been pointed out in Ref. 


24 1 that t 


le bracket adopted in Ref. 


18| is not Poisson, because it 


22,l23j. Very disappointedly, if the Hamiltonian splitting 


does not satisfy the Jacobi identity 

method proposed in Ref. jl^ is applied with the correct Poisson bracket, a.k.a. the Morrison- 


Marsden-Weinstein (MMW) bracket 


22 


25| . one of the subsystems cannot be solved exactly 


in general. Therefore, Hamiltonian integration methods for the Vlasov-Maxwell equations 
cannot be constructed by using the splitting method developed in Ref. j^. 

In this paper, we apply the MMW Poisson bracket 22, l23| and develop a family of 


Hamiltonian integration methods via a new splitting of the Hamiltonian functional. In 
the current context, a Hamiltonian integration method is defined as a numerical integrator 
that preserves the Poisson structure of the Vlasov-Maxwell system specified by the MMW 
bracket. In addition to splitting the Hamiltonian into electrical energy, magnetic energy, and 
kinetic energy, we further split the kinetic energy into three Cartesian components. It turns 
out that exact solutions for all of the resulted subsystems can be calculated by the method 
of characteristics. By compositions of the exact solutions to the subsystems, methods of 
theoretically arbitrarily high order in time for the original Vlasov-Maxwell equations can 
be constructed. The exact solutions of the subsystems preserve the same Poisson structure 
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as the Vlasov-Maxwell equations, so do the combined algorithms. Therefore, the good 
properties of a Hamiltonian integrator, such as the long-term stability and accuracy and 
global bound on the energy error, are all inherited by this family of new integration methods 
for the Vlasov-Maxwell system. 

The Vlasov-Maxwell equations considered in the present study are 


df df df 

_ + ,._ + (E + v><B).- = 0. 


V X B = / v/(x, V, f)dv -h 


aE 

'm' 


V X E = 


OB 


dt ’ 

V • E = /"/(x, v,f)dv, 


V-B = 0, 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


where /(x, v,f) is the distribution function of position x G f/ C and velocity v G at 
time t, and (E(x, t), B(x, f)) G x are the electromagnetic helds. For easy presentation. 


the species index, charge, mass, and other constant are omitted. Equations ([I])-© are c 
and equations (111)-(I5]) result from the gauge symmetry of the system. According to 
VM equations ([I])-© are equivalent to the Hamiltonian system 

dT 


osed. 


23|. the 


dt 


= {{XH}} 


( 6 ) 


on the phase space MV = {(E, B)|V ■ B = 0, V • E = f fdvj, with -F(/, E, B) being any 
functional on MV and l-i being the Hamiltonian functional dehned as 


'H(/,E,B) = ^ J |v|V(x,v,f)dxdv+ ^ J (|E(x,f)|2 + |B(x,t)p) dx. 
Here, {{•,•}} denotes the Morrison-Marsden-Weinstein bracket 


(7) 


21 


123 , 


{{J-,a}}(/,E,B) = j 


Je 


V X 


SB 


S_^ 

(5E 


V X 




dx 


( 8 ) 


(5E 9v Sf (5E dv Sf 


dxdv 


d ST d sg\ . , 
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In the first term on the right-hand side, {•, denotes the canonical Poisson bracket for 
functions of (x, v). With the initial conditions /(x, v, 0) = /o(x, v) and (E(x, 0), B(x, 0)) G 
M.V, there exists unique solution to the system ([6]), on which the bracket ([8]) is Poisson and 
is preserved. Given this Hamiltonian formulation of the Vlasov-Maxwell system, Poisson 
structure preserving integration methods can be constructed as follows. The system is 
hrst split into several solvable subsystems by decomposing the Hamiltonian functional ([7]). 
We then hnd the exact solutions to the subsystems, and hnally the exact solutions of the 
subsystems are composed in a proper way to construct integrators for the Vlasov-Maxwell 
equations that preserve the Poisson structure fISl) . 

Firstly, we follow Ref. fl^ to split the Hamiltonian into three parts. 


l-L = He + Hb + 


r \^) 

/ |v|^/(x, V, f)(ix(iv. 

Using the MMW bracket ([8]) and the Hamiltonian equation ([6]), the VM equations (fT |) -(|5 |) 
can be split into three subsystems on AiV, 


'He = ^ j |E(x,f)|^(ix, nB = ^J |B(x,t)|^dx, Hj 


A={{A,He}}, A = {{A,Hb}}, A = {{A,Hf}}. 


Next, we try to solve the subsystems for exact solutions. The subsystem A = {{-F, 
associated with the Hamiltonian He is equivalent to 


dt 


E(x,t) 

aE _ 


d'v 

0 , 


aB 


-V X E. 


0 , 


( 10 ) 


Given the initial functions /o and (E(x, 0), B(x, 0)) G M.V , the solution to the subsystem 
(UnD is 

/(x,v,t) = /o(x,v-tE(x,0)), 

E(x,t) = E(x,0), (11) 

B(x, t) = B(x, 0) — tv X E(x, 0). 

Using the notation exp(?fEt) as the update operator in time, we denote this solution formally 
as 

(/(x, V, t), E(x, t), B(x, t))^ = exp('HEt)(/o(x, v), E(x, 0), B(x, 0))'^. 
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The subsystem = {{J^, "Hb}} corresponding to the Hamiltonian I-Lb is equivalent to, 

^ = 0 
at ’ 

^ = V X B, (12) 

at 

ae 

ar = °' 

With initial conditions given on WIV, the solution to the subsystem (IT^ is 


= /o(x, v), 

E(x, t) = E(x, 0) + tv X B(x, 0), 
B(x, t) = B(x, 0). 


(13) 


We denote the corresponding update operator as exp{'HBi)- 

For the Hamiltonian "H/, however, the subsystem W = {{W, 'H/-}} is in a more complicated 
form. 


df 9f . / \ \ <9/ 

_ + v.- + (vxB(x,l)).- = 0. 


aE 




(14) 


aB 

'm 


= 0 . 


It has been shown that unless the magnetic field B vanishes or is uniform in space, this 
system can not be solved exactly 2^. Therefore, we search for a Poisson structure preserving 


method for Eq. fflTj) by further splitting the Hamiltonian Hf into more solvable subsystems. 
Utilizing the structure of the cross product in the Cartesian frame. 


V X B = Bv = 


0 i?3 —i?2 

— i?3 0 Bi 

B 2 —Bi 0 


we further split the Hamiltonian into different Cartesian components, i.e., 
W/= Wi/+ W 2 /+ %/, W./= iy t)f/(x,V,()<ix(iv, 1 = 1,2,3. 


(16) 


The subscript of B and v denotes the corresponding Cartesian component of the vector. 
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The subsystem associated with each Hamiltonian 1-iif is 


df ^ , df 

+ - Bi-l{x.)Vi- - 

ot UXi 


Bi+i[yi)vi- -= 0, 

OVi-l 


dEi 

dt 


ni/(x, v,t)dv, 


( 16 ) 


as _ dEi_i _ dEi+i _ 

dt ~ ' dt ~ dt ~ ' 


with ^4 := Ui and Vq := V 3 . The hrst equation of (IT^ can be solved by the method of 
characteristics, and the characteristic equations are 

Xi = Vi, Vi+i = -Bi_i{x.)vi, Vi-i = Bi+i{x.)vi, 


It is a system of ordinary differential equations, whose exact solution is 

Xi{t) = Xi(0) + tvi, 
rxi{t) 

Vi+i{t) = ni+i(0) - / Bi_i{yd)dxi, 

J Xi{0) 
rxi{t) 

Vi-i{t) = ni_i(0) + / Bi+i{x.)dxi. 

J Xi(0) 

Therefore, for the initial conditions /(x, v,0) = /o(x, v) and (E(x, 0), B(x, 0)) G A4.V, the 
exact solution to the subsystem (1T6|) is 

/(x, V, t) = /o(x - tViGi, V + Fi_iei +1 - Fi+iei_i), 


F, = 


Bi{'id)dxi, / = z + 1, z - 1, 


' Xi—tVi 



(17) 


E(x,f) = E(x, 0) - I I eiVif{y.,\,T)dYdT, 


B(x,f) = B(x,0),, 

where is the unit vector in the Tth Cartesian direction. We denote the exact solution 
given by Eq. (fTT)) as (/, E, B)(f)^ = exp('Hi/t)(/o(x, v), E(x, 0), B(x, 0))^. More specihcally, 
the exact solution of the subsystem corresponding to T-Lif is 

/(x, V, t) = /o(xi - tvi, X2, X3, Vi, V2 + V1E3, V3 - V1F2), 

Fi = f Bi{x.)dxi,l = 2,3, 

J xi—tvi 

Ei(x,t) = Ei(x,0) - ^ J ni/(x, v,r)dvdT. 
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The exact solution associated with exp(7^2/0 is 


/(x,V,t) = fo{Xi,X2 -tV2,X3,Vi - V2F3, V2, V3 + V2F1), 
rX 2 

Fi= Bi{x.)dx2,l = 1,3, 

J X2—tV2 

^ 2 (x,t) = ^2(x,0) - ^ j U 2 /(x, V,T)dvdT, 


and the exact solution associated with exp{'H 3 ft) is 


/(x, v,t) = fo{xi,X2,X3 - tV3,Vi + ^3^2, ^2 - n3Fi,n3), 
r^3 

Fi= Bi{x)dx3,l = 1 , 2 , 

J X3—tV3 

^3(x,t) = E3(x,0) - ^ J n3/(x, v,T)dvdT. 


Overall, we split the Vlasov-Maxwell equations into five subsystems. 


di — Be + Bb + Bif + B2f + B3f, 

F = {{F,Be}}, F = {{F,Bb}}, F={{F,Bif}}, F = {{F,B2f}}, F={{F,B3f}}. 

with Be,Bb and Bif, i = 1,2,3, defined in Eqs. ([9]) and (IT^ . and the Poisson bracket 
being defined in Eq. ([H]). The exact solutions to the subsystems are given explicitly by 
Eqs. flTT]) . flT^ . and flT7|) . respectively. For initial functions defined on MV, the solutions of 
the subsystems are all on AIV and preserve the MMW Poisson bracket. 

Given exact solutions to the subsystems, integration methods for the original Vlasov- 
Maxwell equations can be constructed by compositions. For example, if we denote the 
solution at time t as Z{t) := (/(x, v, t), E(x, t), B(x, t))^, a first order numerical update in 
one time step with step size At can be derived from the Lie-Trotter composition 


Z{t + At) = exp{AtBE) exp{AtBB) exp{AtBif) exp{AtB 2 f) exp{AtB 3 f)Z{t), 

and a second order symmetric method can be constructed by the following symmetric com¬ 
position, 

Z{t + At) =exp exp ^xp ^xp exp (AtBsf) 

exp exp exp exp Z{t). 
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With the help of the BCH formula, proper compositions can be found to yield methods 


of arbitrarily high order in time 


26| . Given the geometrical properties of the solutions 


to the subsystems, the combined methods will preserve the Poisson structure and generate 
solutions on M.V, if the initial values Z{0) are taken on M.V. 


In conclusion, we have constructed a family of Poisson structure preserving, or Hamilto¬ 
nian, integration methods for the Vlasov-Maxwell equations by the splitting technique. The 
Hamiltonian is split into hve parts, each part associates with a solvable Hamiltonian sub¬ 
system with respect to the MMW Poisson bracket. As a consequence, integration methods 
for the Vlasov-Maxwell equations constructed via composition of the exact solutions of the 
subsystems will preserve the original Poisson structure. These new Hamiltonian methods 
are expected to exhibit long-term accuracy and hdelity, as well as bounded error in energy 
and other invariants. Numerical applications of the methods will be reported in future 
publications. 


In the present study, we have focused on the time integration of the Vlasov-Maxwell 
system and the preservation of the original Poisson structure. If the system is discretized 
in space, the resulted scheme needs to preserve a discrete Poisson bracket depending on the 
spatial discretization scheme 19| . 


As a hnal note, we emphasize that this general methodology of constructing Hamiltonian 
integration algorithms via Hamiltonian splitting is also applicable to other systems that 
admit Poisson structures. Recently, Burby et al. discovered the Poisson structures for the 
gyrokinetic system j^| and the collision operator j^|. Application of the splitting technique 
to these systems will generate more effective numerical algorithms for large-scale, long-term 
simulation studies of plasma physics. 
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